Magnetism in parent Fe-chalcogenides: quantum fluctuations select a plaquette order 
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We analyze magnetic order in iron-chalcogenide Fei+yTe - the parent compound of high- 
temperature superconductor Fei+yTei_a;Sea;. Neutron scattering experiments show that magnetic 
order in this material contains components with momentum Qi = (7r/2, tt/2) and Q2 = (7r/2, —it/2) 
in Fe-only Brillouin zone. The actual spin order depends on the interplay between these two com- 
ponents. Previous works argued that spin order is a single-Q state (either Qi or Q2). Such an 
order breaks rotational C4 symmetry and order spins into a double diagonal stripe. We show that 
quantum fluctuations actually select another order - a double Q plaquette state with equal weight 
of Qi and Q2 components, which preserves C4 symmetry but breaks Z4 translational symmetry. We 
argue that the plaquette state is consistent with recent neutron scattering experiments on Fei+yTe. 



Introduction. The analysis of magnetism in parent 
compounds of iron-based superconductors (FeSCs) is an 
integral part of the program to understand the origin of 
superconductivity in these materials . Parent com- 

pounds of Fe-pnictides are moderately correlated met- 
als, whose resistivity increases with increasing T, and 
the electronic structure is at least qualitatively consis- 
tent with that of free electrons on a lattice 0, • Mag- 
netic order in such systems can be reasonably well under- 
stood within itinerant scenario 0-0, Q due to enhance- 
ment of free-electron susceptibility at momenta connect- 
ing hole and electron Fermi surfaces (FSs). The loca- 
tions of the FSs select two possible momenta for the 
order - (0,7r) and (7r,0)- in the Fe-only Brillouin zone 
(BZ). Electron-electron interaction and the shape of the 
FSs further reduce the ground state manifold to single- 
momentum states with either (0,7r) or (7r,0), but not 
their mixture [§] ) . In each of these two states spins are 
ordered in a stripe fashion ~ ferromagnetically along one 
direction in 2D Fe-plane and antiferromagnetically in the 
other. Such an order breaks C4 lattice rotational sym- 
metry and causes pre-emptive spin-nematic order jlSl ). 
The same magnetic order is selected in the strong cou- 
pling approach, which assumes that the system is not 
far from Mott transition, and magnetic properties are 
reasonably well described by Ji — J2 model with near- 
est and second-nearest neighbor spin exchange 
The actual coupling in Fe-pnicties is neither truly small 
nor strong enough to cause Mott insulating behavior [l^ , 
which makes it extremely useful that the two descriptions 
agree. Upon doping, long-range order is lost, but mag- 
netic fluctuations evolve smoothly and remain peaked at 
or near (0, tt) and (tt, 0) even beyond optimal doping [l8l |. 

There is one family of FeSCs - 11 Fe-chalcogenidcs 
Fei+yTei-ajSe^;, in which smooth evolution between par- 
ent and optimally doped compounds does not hold. Mag- 
netism in these materials changes considerably between 
X = and x ~ 0.5, where the Tc is the largest. Near 
optimal doping magnetic fluctuations are peaked at or 
near (0,7r) and (7r,0), as in Fe-pnictides, while mag- 



FIG. 1: The two possible coUinear configurations for the 
Ji — J2 — J3 model: (a) orthogonal double stripe (ODS) and 
(b) diagonal double stripe (DDS). 



netic order in a parent compound Fei+yTe has very dif- 
ferent momenta ±(7r/2, ±7r/2) |19l - [23| . Upon doping, 
the spectral weight at ±(7r/2, ±7r/2) decreases, and the 
spectral weight at (0,7r) and (7r,0) increases [20[. The 
transport properties of Fci+^Te are also quite different 
from those of parent compounds of Fe-pnictides: the 
resistivity, p(T), of Fci+^Te does not show a promi- 
nent increase with increasing T, but instead remains flat 



and even shows a small increase as T decreases [24 1. 
ARPES studies of Fci+^Te show that low-energy spec- 
tra are very broad psj . consistent with the notion that 
electrons are not propagating. These observations lead 
several groups to suggest that parent Fe-chalcogenides 
are more correlated than parent Fe-pnictides, and mag- 
netism in Fei_|_j,Te can be understood by assuming that 
electrons are "almost" localized and interact magneti- 
cally via a Heisenberg exchange 26|-29|. This scenario 
is in line with a more generic idea |30l - l32| that in any 
FeSc, a certain percentage of electronic states are local- 
ized and phase separated from itinerant electrons, and 
the percentage of localized states varies between differ- 
ent materials. An alternative scenario for FeTe, which 
we don't discuss here, is orbital order [s^l 

In this communication we apply the localized elec- 
tron scenario to Fei+yTe and verify whether the observed 
commensurate ±(7r/2, ±7r/2) order can be obtained in a 
Heisenberg model with exchange interactions up to third 
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neighbors. Classically, ±(7r/2, ±7r/2) order is unstable 
with respect to a spiral order for any non-zero first neigh- 
bor exchange, unless one artificially breaks C4 symm etry 
and sets interactions to be spatially anisotropic 2l|, |30| . 
We analyze the isotropic quantum Heisenberg model and 
show that quantum fluctuations do stabilize a commen- 
surate ±(7r/2, ±7r/2) order in some range of parame- 
ters. However, this stabilization does not uniquely de- 
termine spin configuration as a generic ±(7r/2, ±7r/2) 
order is a superposition of two different Q— vectors: 
Qi = (7r/2,-7r/2), and Q2 = (7r/2,7r/2): S(r) ^ 
Ai cosQir-f- A;^ sinQir + A2 cos Q2r-|- Aj sinQ2r, with 
|A,| |A,^| = A and Aj • A2 = A[ ■ A'^ = 0. In Fig. [T] 
we show two prototypical commensurate spin configura- 
tions - a single Q bi-coUinear spin order (Ai = A^ = A, 
A2 = A2 = 0), which breaks C4, and a double Q pla- 
quette order (A'^^ = A2 = A, Ai = Aj = 0), which 
preserves C4 symmetry, but breaks Z4 translational sym- 
metry (four equivalent plaquette states are obtained by 
moving a black square in Fig. [TJi by one lattice site in 
either x or y direction). Bi-collinear spin order is often 
called diagonal double stripe (DDS), and plaquette or- 
der is called orthogonal double stripe (ODS), and we use 
these notations below. The real-space configuration for 
both orders is " up- up-down-down" along x and y direc- 
tions. 

Most of previous theoretical and experimental works 
assumed that the commensurate order is DDS [l^l and 
studied in detail the feedback from this order on elec- 
We argue that quantum fluctuations of spins 



21 



trons 

actually select ODS order as a stable collinear state for 
weak but finite nearest- neighbor exchange Ji , while DDS 
state is unstable for any non-zero Ji. The DDS and 
the ODS orders have qualitatively different forms of the 
static structure factor S{q) (two peaks vs four peaks), 
but this is difficult to detect in real materials because 
of domains. The authors of [i^l however argued that 
form of S{q) in a paramagnetic phase allows one to dis- 
tinguish between DDS and ODS, even in the presence of 
domains, and found that their results are consistent with 
strong ODS fluctuations. Another argument in favor of 
the C4 preserving ODS spin order is the absence of or- 
thogonal distortion in Fci+yTc - there is a monoclinic 
distortion below Tn, but this docs not break rotational 
in-plane C4 symmetry. There is also numerical evidence - 
ODS order has been found in exact diagonalization stud- 
ies of Jl — J2 — J3 model on clusters up to 36 spins [s^ . 
The same ODS order has been found in the mean-field 
studies of the t — J model in another Fe-chalcogenide 
K0.8Fe1.6Se2 M- 

Model. We follow earlier works and model magnetic 
interactions in Fei4_JTe by a Ji — J2 — J3 Heisenberg 
model [23,[2i[3lii: 




H 



11=1 {ij) 



(1) 



FIG. 2: Spin order in the classical Ji — J2 — Js model at Ji = 
0. Classically degenerate configurations form four sublattices, 
labeled as a, b, c, and d. A configuration with arbitrary 7, 6, 
and is a ground state. In our notations, sublattice spins are 
Ai + A2, Ai - A2, A'i — A2, and A^ + A2, respectively. 



where Ji , J2 , and J3 are antiferromagnetic exchange cou- 
plings between first-, second-, and third- nearest neigh- 
bors. For Fei+yTe the values of Ji,J2, and J3 have 
been estimated in (27| and found to be in the range 
J3 > ^ » Jl. In this limit, the classical ground state 
of ([1]) is a spiral with the pitch vector Q = (±g, ±g), 
where q = arccos( jj~^;^j^ ) [37|. At Ji = 0, the model 
has an extensive degeneracy, and any order with momen- 
tum ±(7r/2, ±7r/2) is the classical ground state, includ- 
ing DDS, ODS, and an infinite number of other four- 
sublattice states (Fig. 

We consider here what happens in the quantum model, 
at a finite Ji. We show that the ODS state is unambigu- 
ously selected by quantum fiuctuations to be the ground 
state in some range of Ji, before a spiral order sets in. 
Our key reasoning is that only some classically degener- 
ate ground states at Ji = are degenerate by symmetry; 
others are "accidentally degenerate". The situation is 
quite similar to the one in the well-known Ji — J2 model 
at J2 > Ji/2 [3I]. We argue that quantum fluctuations 
lift accidental degeneracies and gap out some of the spin- 
wave modes which in the classical limit become unstable 
(imaginary) at Ji 7^ 0. For the DDS state the lifting of 
the accidental degeneracies does not help, as the modes 
which become unstable at a flnite Ji are the true Gold- 
stone modes at Ji = 0. On the other hand, for ODS state 
classically unstable modes are accidental zero modes at 
Jl = 0, and quantum fluctuations lift the energies of 
these modes to flnite values, making ODS the state sta- 
ble in a flnite range of Ji. We verified that ODS state is 
indeed the ground state in this range. 

Large-S spin-wave calculations. We consider large 
value of spin S and study the role of quantum fluctu- 
ations within 1/5 expansion. The computational steps 
are presented in [s^. For J^ = 0, spins on even and odd 
sites form two non-interacting sublattices, each described 
by J2 — J3 model. This model is identical to " Ji — J2" 
model, with diagonal hopping J2 playing the role of " Ji" 
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and third-neighbor hopping J3 playing the role of " J2". 
One can use this analogy and borrow the results of the 
quantum analysis of " Ji — J2" model [s^l- For J3 > J2/2 
(which holds in Fei+^Te), quantum fluctuations select 
stripe configurations within each sublattice, i.e. the an- 
gle 7 in Figl2]is locked at 7 = or 7 = tt, and the angle 
9 is locked &% 9 = (f) ot: 9 = -k + (p. The states with 
7 = and 7 = tt are equivalent up to an interchange 
of X and Y directions, and below we set 7 = 0. The 
coUinear DDS and ODS states belong to the manifold of 
selected states and correspond to different locking of the 
angle <j) between the nearest-neighbor spins: DDS state 
corresponds to i/i = 0, ^^ = tt or = tt, = 0, while ODS 
corresponds io(t) — 6 — QoT:(f) = 9 = Ti. 

To analyze whether a generic state selected by quan- 
tum fluctuations at Ji = remains stable at a finite 
value of Ji , we need to know its excitation spectrum. At 
Ji = 0, spins on even and odd sites are decoupled, each 
sublattice is described by its own bose field (ak for even 
sites and /3k for odd sites) , and spin- wave excitations are 
described by 

= 5(l]„k4ak + f^^ik/3i^/3k), (2) 

The classical spin-wave spectrum is the same for all se- 
lected states 

fik = S{Al - Blf/^, Ak = 4J3 + 2J2 cos(fc, + ky), 
Bk = 2J2(cos2fc2; -f cos2fcj,) +2J2Cos(fc^ — ky). (3) 

This spectrum contains nodes at ±(7r/2, ±7r/2), but some 
of them are not symmetry-related and are lifted by quan- 
tum fluctuations. For the sublattice made of even sites, 
the order has momentum ±(7r/2, — 7r/2) (Fig. [2]b), hence 
the true nodes are located only at these momenta, while 
the ones at ±(7r/2,7r/2) must be lifted. For the sub- 
lattice made out of spins at odd sites, the order has 
momentum ±(7r/2,7r/2) if we take = 0, like in the 
ODS, and ±(7r/2, -7r/2) if we take 6* = tt + 0, like in 
the DDS. Quantum fluctuations then must lift the nodes 
at ±(7r/2,-7r/2) and at ±(7r/2,7r/2) for the ODS and 
the DDS state, respectively. We computed quantum cor- 
rections to the spectrum in Eq. ^ within perturbation 
theory to order 1/S and indeed found that accidental 
nodes are lifted by quantum fluctuations and only true 
Goldstone modes remain [soj . 

We next set Ji to be small but finite and consider 
which of stripe states, if any, remain stable. The quali- 
tative reasoning is the following: a non-zero Ji couples 
the two sublattices and adds to the Hamiltonian ([2]) the 
terms in the form aj^/3k and ak/3k- For the DDS state (or, 
more accurately, for the DDS family of states as we keep 
(/) as a parameter) the stripes on even and odd sites are 
directed parallel to each other, and the dispersions of ak 
and /3k fields are identical, including 0{1/S) terms. The 
two dispersions are then gapless at the same momenta 



k = ±(7r/2, — 7r/2). Around these k points, the pertur- 
bation theory in Ji is singular, as there is no symmetry 
requirement which would force the coupling to vanish 
at ±(7r/2, — 7r/2). As a result, the excitations become 
purely imaginary close enough to ±(7r/2, — 7r/2), which 
implies that the DDS states are unstable at any non-zero 
Ji. On the other hand, for the ODS family of states, 
the dispersions Q,^ and Q,^ have nodes at different mo- 
menta, ±(7r/2, — 7r/2) and ±(7r/2, 7r/2), respectively. Be- 
cause of this disparity, perturbation theory near either 
±(7r/2, — 7r/2) or ±(7r/2,7r/2) is not singular, and correc- 
tions in Ji only gradually shift the values of spin-wave 
velocities thus keeping ODS states stable. 

We verified this reasoning by explicit calculations. We 
first obtained the Ji-induced interaction in terms of the 
original Holstein-Primakoff bosons and then re-expressed 
it in terms of ak and /3k bosons from Eq. ([2]) , which are 
related to the original ones by Bogoliubov transforma- 
tion. The UkWk-coefHcients of this transformation dress 
up the interaction terms. For the DDS states, expand- 
ing the Hamltonian near the true Goldstone points at 
(7r/2, — 7r/2) as k = (7r/2, — 7r/2) -I- k we obtain Hdds = 
Hsiv + SHdds, where Hsw is given by ([2|) with 

ni^n^^ 4S^M2J3 + J2)Ckl + kl- 2aWkyf'\ (4) 

where a = < 1, and 

6Hdds = A^^'5k(al/3£ + a^^^^^ + h.c) (5) 

where 

^DDS 2^3 + ~ky - h .g^ 

^ 2 \ J3 J (kl^kl-2ahkyYn 

The coupling term remains finite when k^^y tends to zero, 
except for special directions. Diagonalizing ^ we find 
that at low enough k one of the two solutions is i?? k, 

-2r2?^'^'A?^'S'. A negative -B? implies that fluctuations 
around a DDS state grow exponentially with time and 
make this family of states unstable. 

For the ODS states the situation is different because 
near any of the points ±(7r/2, —tt/2) or ±(7r/2, 7r/2), the 
zero in one of the spin- wave branches is lifted by quantum 
fluctuations. For example, near (— 7r/2,7r/2) expanding 
of the Hamiltonian again gives Hqds — Hsw + SHqds, 
however now only fi^ is gapless, while il^ is gapped with 
the gap of the order 1/S. The interaction term 5 Hqds 
has the same form as in but with 

A?^^ = 2Ji52(2j ^ o ( Ji53/2|fc|i/2A 

^ i V ^ ^(r2?r2^)i/2 V ' ' y 

(7) 

Diagonalizing Hqds we find two solutions, 

EI2 = l{inir + inir (8) 
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One of the solutions is gapped to order 1/ the other is 
linear in k with a stiffness which differs from its value at 
Ji = by O^JiS/Ja). We see that the ODS states are 
stable (for any (j)) as long as J1S/J3 is small. 

On a more careful look, we find that the ODS spin or- 
der allows for Ji— induced umklapp processes, which also 
renormalize the dispersions of the ODS states. Indeed, 
because ODS state breaks Z4 translational symmetry, 
the Ji interaction contains not only the terms at zero 
transferred momentum, as in ([5]), but also terms with 
momentum transfer in multiples of tt along each axis. 
Near k = (7r/2, — 7r/2), the most relevant of such umk- 
lapp terms is the one with momentum Q = (0, tt), which 
connects a gapless boson at {tt/2, —tt/2), and a gap- 
less boson at (7r/2,7r/2). However, because breaking 
of Z4 is equivalent to breaking local inversion symmetry 
(a reflection around one column or one row in Fig. H^), 
the umklapp vertices /^^''-'^^ contain extra momentum 
gradient compared to non-umklapp vertices. In explicit 
form, we find at small k = k — (7r/2, —tt/2), 



the C4 symmetry of the static structure factor S{q) ob- 
tained in neutron scattering experiments (22j . We com- 
puted S{q) for both the DDS and the ODS states, and we 
indeed found that the structure factor for the ODS or- 
der has four identical peaks at (±7r/2, ±7r/2), while the 
structure factor for the DDS state has only two peaks at 
{tt/2, -tt/2) and {-tt/2, tt/2). While the observed four 
peaks are consistent with ODS, we caution that the ab- 
sence of the anisotropy in the structure factor obtained 
in neutron scattering could be due to the twinning of the 
crystal. However, as the magnetic domain's structure of 
the crystal can be controlled using polarized neutrons, 
the careful analysis of the neutron scattering data might 
dissect the contribution from different domains. The au- 
thors of Ref. [2^ made another argument that, even in 
a twinned crystal, the form of S{q) throughout the Bril- 
louin zone differentiates between strong DDS and ODS 
fiuctuations, and argued that their data are more con- 
sistent with tendency towards ODS order. This again 
agrees with our results. 



^U,ODS 



.Ji 



1/2 



k+Q 



(9) 



where f]?,!!? 



45(J3(2J3 ± J2))'/Hkl 



2kxky)^/'^ and the angle (j) specifies the spin order within 
the ODS family of states. We see that /s^'-'^^ scales lin- 
early with k, i.e., is of the same order as ^'jf ■ We com- 
puted the corrections to spin-wave velocity and found 
that they scale as Ji/a/4J| — Ji, i.e., are small. At 
the same time, we see from the Eq.(|9]) that de- 
pends on the angle (j). Respectively, the corrections to 
the ground state energy also depend on and should se- 
lect which state within the ODS family has the lowest 
energy. The computation is straightforward and yields 
AEgr — —A cos^ (p, with A > 0. We see that the collinear 
ODS state, for which = or tt. is indeed the state with 
the lowest energy. 

The outcome of our analysis is that the collinear ODS 
state remains stable and has the lowest energy within a 
family of similar states. At small Ji, the ODS state has 
a finite stiffness towards fiuctuations which tend to break 
collinear order in favor of a spiral one. The ODS state 
remains stable up to Ji ~ J3/S, at larger Ji the stiffness 
changes sign, and the system develops a spiral order. 

Experimental signatures of ODS state. Because the 
ODS state does not break C4 translational symmetry, it 
does not cause a pre-emptive spin-nematic order, in con- 
trast to parent compounds of other FeSCs The data 
for Fei+yTe show that the system develops a monoclinic 
distortion below a certain T, but in-plane C4 symmetry 
remains unbroken (it only breaks in doped compounds 
Fei+yTei_a;Sea; with x > 0.5 0). The unbroken C4 
symmetry in the ordered state also manifests itself in 



Summary. In this communication we analyzed the 
type of magnetic order in Fei+y Te - the parent compound 
in a family of Fe-chalcogenide superconductors. The 
magnetic order in this material is different from in other 
parent compounds of FeSCs - spins are ordered in up- up- 
down-down fashion (Fig. [5]). Experiments show [2J, |40[ 
that the tendency towards Mott physics is stronger in 
Fei+j,Te than in other parent compounds of FeSCs, sug- 
gesting that the magnetic order in Fei+j,Te can be rea- 
sonably well understood within the localized scenario by 
solving the Heisenberg model with exchange interaction 
extendin g up to third neighbors [27[. Several groups ar- 
gued [Hlill, 113, m that the ordered up- up- down- down 
spin configuration is diagonal double stripe. Such an or- 
der breaks C4 lattice rotational symmetry. We argued, 
based on our analysis of quantum fiuctuations in the 
Heisenberg model with first, second, and third- neighbor 
exchange, that such a state is unstable, but another up- 
up-down-down state - the orthogonal double stripe, is 
stable and is the ground state in some parameter range. 
This state (which is also called a plaquette state) breaks 
Z4 translational symmetry, preserves C4 symmetry, and 
does not cause orthorhombic distortion. Also, its struc- 
ture factor has four equivalent peaks at (±7r/2, ±7r/2), 
in agreement with recent neutron scattering studies of 
Fei+yTe An interesting issue that deserves further 
study is whether Z4 translational symmetry can be bro- 
ken before a true ODS spin order sets in, as it happens 
in other systems 42 1 . 
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In this Supplementary Material we provide details of 
our analysis of how quantum fluctuations lift accidental 
nodes in the spin-wave spectra of the sub-set of states 
selected by quantum fluctuations in the Ji — J2 — J3 model 
at J3 > J2/2 and Ji = (i.e., in the J2 — J3 model with 
^3 > <^2/2). The Hamiltonian for the model is given by 
Eq. (1) in the main text. At Ji = spins on even and 
odd sites do not interact with each other. It is therefore 
sufficient to consider only the sub-set made out of spins 
on, say, even sites. 

Classically degenerate ground states are shown in Fig. 
2 of the main text. For spins on even sites, any two- 
sublattice configuration with arbitrary angle 7 between 
spins in a and b sublattices is the ground state. Like we 
said in the text, quantum fluctuations break this degen- 
eracy and select the states at 7 = or 7 = tt as the only 
two ground states. The state with 7 = can be viewed 
as a spiral state with Qi = (7r/2, — 7r/2), and the 7 = tt 
state corresponds to a spiral with Q2 = {it/2,tt/2). Be- 
cause only even cites are involved, configurations with Qi 
and —Qi are identical. 

The classical spectrum of either Qi or Q2 state con- 
tains zero modes at corresponding Q, which must be 
there by Goldstone theorem, but also contains zero 
modes at "wrong Q" (at Q2 for the Qi spiral and vise 
versa) . These last modes are not associated with symme- 
try breaking (i.e., are "accidental") and must be lifted by 
quantum fluctuations. In the main text we stated that 
quantum fluctuations do gap accidental zero modes and 
studied the consequences. Here we show how this actu- 
ally happens. For definiteness, we consider Qi configu- 
ration (7 = 0). 

The lifting of accidental zero modes can be studied 
within the original two-sublattice picture, with two dif- 
ferent bosons describing fluctuations of spins in a and b 
sublattices. It is more convenient, however, to perform a 
uniform rotation to a local reference frame in which the 
magnetic order is ferromagnetic and describe the exci- 
tation spectrum using just one Holstein-Primakoff bose 
operator a^: 



Sf = 5 - aja, 
S+ = {2S-ala,)ia, 
5r = al(25 



(S.l) 



As it is customary to spin- wave analysis, we assume 
that S is large. Our goal is to obtain the excitation spec- 
trum for Qi configuration including terms of order 1/5*, 
which describe quantum fluctuations. Accordingly, we 



substitute Eq. (S.I) into Hamiltonian (Eq. (I)) and then 
expand it to the quartic order in a field. We obtain 



H = Eo + S{H2 + Hi) , 



{S.2) 



where H2 and Hi are the quadratic and the quartic 
terms, respectively. In momentum space, we have 

H2= Ek(^k4«k + f^(aka-k + /i.c.)), (5.3) 

where 

= U:i + 2J2Cos{k^ + ky) [SA) 
Bk = — 2 J3(cos(2fc^) + cos(2fcj,)) — 2 J2 cos{kx — ky) 

and k is defined in the first magnetic BZ. The quartic 
part of the large- S expansion is given by 



Hi 



J2 



8NS 



^^(^kjj aLi^^ak2ak3ak4 cos(fc4:j - fe 



-alki'^-ks'^kaaki cos(fc4j; 4- kiy) (5.5) 
-2al^^a^^^ak3ak4 sin(fc2a- -I- ki^) sm{k2y + kiy) + h.c. 



jf^kjj aLkiak2ak3ak4(cos2fc4a; 4- cos 2A:4.y) 



8NS \^ 

-a^^.^a^^^ak3aki(cos2(fc2y -I- kiy) + cos2(/s2k -I- ki^)) + h.c. 

The quadratic Hamiltonian H2 is diagonalized by in- 
troducing the Bogoliubov transformation, ctk = WkOk — 
''^kQ'LkJ where the coherence factors Uk and v-^ arc deter- 
mined by 



Uk 



Wk = 



1 / ^k + »k 

2 V nv 



signSk / A k - ^^k 



and 



The diagonalized Hamiltonian H2 is given by 
H2 = E2 + Y, r!ka],ak , (5.6) 



where 
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E2 = ^Y.^-A^ + n^) 



(5.7) 




FIG. 3: Excitation spectrum for Qi spiral. Top panel: The 
spectrum of non-interacting bosons, Eq.(S.6) This spectrum 
contains zero modes at ±(7r/2,7r/2) and ±(7r/2, — 7r/2) Bot- 
tom panel: the spectrum renormalized by quantum fluctua- 
tions to order 1/S. The true Goldstone modes at ±Qi = 
-'-(§' ~f) remain, but the accidental zeroes at ±Q2 = 
±(^, ^) are removed by quantum fluctuations. 



is the contribution to the ground state energy from non- 
interacting a-bosons 

The spin-wave dispersion of non-interacting bosons is 
f^k- At k = ±Qi and k ±Q2, Aq^^^^ = Bq^^^^ = 
4J3±2J2, where -f- and - signs are for Qi and Q2, respec- 
tively. As a result, fig^ = fig^ = 0. However, as we said, 
only the zero mode at k = ±Qi is the true Goldstone 
mode, the other one is accidental and must be gapped 
by quantum fluctuations. To show this, we have to com- 
pute 1/S corrections to the spectrum at these points. 

The 1/S contribution to fik is obtained by evaluat- 
ing the first-order correction from Hj^ what amounts to 
decoupling of H^^ into the product of two pairs of a— 
bosons and replacing one pair by its average value for 
non-interacting bosons. We first define the averages 

- «aq) = vl , (5.8) 

iq = (aqa_q) = (aLqa]j) = UqUq . 

The decoupled is 

k,q 



/ , t t^-^2(k,q) 
(okO-k + a_kO 7, 



(5.9) 
where 

i^i(k,q) = -2J3(^4A'q(cos2(fc^ - q^) + cos^ {ky - qy 

— Lq(2(cos 2qx + cos 2qy) + cos 2ka: + cos 2ky)j 

-2J2(^Kq{4:aiii{kx - qx) sin(fcj^ - qy) 
+2 cos{kx + ky) +2 cos{qx + %)) 
-Lq{2 cos{qx - qy) + cos{kx - ky))^ , 

F2(k, q) = 2 J3^A'q(cos 2(7i, + cos2qy 

-|-2(cos2fca; -I- cos2fcj^)) 

-2Lq(cos2(fc2; - qx) + cos2{ky - qy))^ 

+2J2(^Kq{cos{qx - qy) + 2cos(fc^ - ky)) 
+Lq(-4sin(fca; - qx) sin{ky - qy) 
- cos{kx + ky) - cos{qx + qy))^ 

Combining the decoupled with H2 expressed 
terms of original a-operators we obtain 



H2 + H4=Y1 [(^i< 
k 



A. 



(4) 

k \ t 

^aj^Ok 



5 



1 



(4) 



where 



{Bk + -^){a^a^k + h.c.) , (5.10) 



= 



At k = Qi and k = Q2 wc have 



A 



(4) 

Ql(2) 



-1 



— J2 J3(4/^q(sin2(g,) + Sin^iqy)) 



(4) 

Qi( 



— Lq(2(cos2(7j; + cos2g.y) — 2)^ 

+J2 {k^{-=\^A cos{qx) cos{qy) 
±2 -I- 2 cos((3'j; + qy)) 
-Lq{2cos{qx - qy) T 2)), 

q 

-|-2Lq(cos 2qx + cos 2qy) 

+J2(Aq(cos((7x - qy) T 2) 
-|-Lq(±4cos(qi.) cos(g.y) 
=Fl - cos{qx + qy)) 
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where the upper and the lower signs correspond to Qi 
and Q2 vectors, respectively. One can easily make sure 
that = Sg], but a'^^I ^ B^^l Thus, when 1/S 

corrections are included, the spectrum at Qi remains 
gapless, as it should be, but at Q2 the gap opens up, i.e., 
quantum fluctuations lift the accidental degeneracy at 
the "wrong Q" . We computed the renormalized spectrum 



for all k and show the results in Fig. 3. The renormalized 
spectrum, shown in the bottom panel, clearly has a gap 
at ±(7r/2, 7r/2), where the spectrum of non-interacting 
bosons (top panel) has zero modes. We verified that the 
same effect (lifting of accidental nodes) holds if we use 
biquadratic spin interaction instead of quantum fluctua- 
tions. 



